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A new interpretation of the optical and energy-loss spectra of hexagonal boron nitride is pro- 
vided based on first-principle calculations. We show that both spectra cannot be explained by 
independent-particle transitions but are strongly dominated by excitonic effects. The lowest direct 
and indirect gaps are much larger than previously reported. The direct gap amounts to 6.8 eV. 
The first absorption peak at 6.f eV is due to an exciton with a binding energy of 0.7 eV. We show 
that this strongly bound Frenkel exciton is also responsible for the low frequency shoulder of the it 
plasmon in the energy-loss function. Implications for nanotube studies are discussed. 

PACS numbers: 78.20.-e, 71.35.Cc, 71.45.Gm 



Hexagonal boron nitride (hBN) is isoelectronic to 
graphite and has a similar layered structure. Recently, 
the electronic structure and the optical properties of 
hBN have regained interest due to the discovery of BN- 
nanotubes 1 ^ which can be considered as cylinders formed 
by rolling a single sheet of hBN onto itself. While 
graphite is a semi-metal with a linear crossing of the tt 
and 7r* bands at the K-point, in hBN the different elec- 
tronegativities of boron and nitrogen lift the degeneracy 
at the K-point and lead to a large gap. It was predicted 3 - 
that hBN has an indirect quasi-particle gap of 5.4 eV and 
a minimum direct quasi-particle gap of 6.2 eV. Measure- 
ments of the optical absorption coefficient^*SiL& yield 
unanimously an absorption peak with a maximum be- 
tween 6.1 and 6.2 eV. The spectra in these papers are 
interpreted in terms of independent electron transitions 
from the valence to the conduction band. The onset of 
the absorption peak (varying in the different experiments 
between 5.2 and 5.5 eV) is consequently considered to be 
a measure of the band gap. However, calculations^ at the 
Random Phase Approximation (RPA) level, i.e., in the 
picture of independent electron excitations, yield spectra 
whose shape is not in good agreement with the experi- 
mental spectra. Furthermore, in a recent luminescence 
experiment 10 , evidence for a direct bandgap at 5.97 eV 
and a series of bound excitons with a maximum binding 
energy of 0.149 eV was found. In this communication, we 
reinterpret the different experimental results and demon- 
strate the importance of excitonic effects for a proper un- 
derstanding of both absorption and electron-loss spectra. 

For the interpretation of the spectra, we calculate the 
real and imaginary parts, e±(E) and 62(E), of the energy- 
dependent macroscopic dielectric function on two levels 
of approximation: i.) using RPA, ii) on the GW+BS 
level, i.e., including quasi-particle corrections on the 
level of the GW-method and calculating excitonic effects 



through the Bethe-Salpeter (BS) equationAi. In both 
cases, we start with a calculation of the wavefunctions 
using density- functional theory (DFT). With the code 
ABINITpiS we calculate 1 ^ the valence-band wavefunctions 
and a large number of conduction-band states, if) n ki with 
a band-index n and a sufficiently fine sampling of the 
crystal momentum k. In this work we use two descrip- 
tions of exchange-correlation effects: the customary lo- 
cal density approximation (LDA) and the exact-exchange 
(EXX) potential. To obtain the EXX potential we fol- 
lowed Ref . O 

On the level of the RPA (which corresponds to the 
use of Fermi's golden rule for 62(E)), we compute the 
optical spectra in terms of the dipole-matrix averaged 
joint-density-of-states. However, ground-state DFT pro- 
vides only a qualitative and, often incorrect, description 
of the excited states. The reason is that many-body ef- 
fects beyond the independent particle approximation are 
neglected. But in systems with a poor electronic screen- 
ing, the electron-hole attraction can be far from negli- 
gible. Good agreement with experiments is restored by 
using the self-energy approach of many body perturba- 
tion theorjiii. Starting from the DFT energies and wave- 
functions, we calculate the quasi-particle energies ( "true" 
single-particle excitation energies), E n k, by solving the 
Dyson equation: 



V_2 

2 



-— + Vext + VHartree + 



V-nk ~~ ^nkVnk- 



(o.i) 

The self-energy E is approximated as the energy con- 
volution of the one-particle Green's function G and the 
(RPA) dynamically screened Coulomb interaction W, 
(E = iGW). We perform a "semi-self consistent" (GW ) 
calculation by updating the quasi-particle energies in G 
(but not in W) until the resulting quasi-particle ener- 
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FIG. 1: Bandstructure of hBN: a) LDA (dotted lines), LDA-GW (circles); b) EXX (dotted lines), EXX-GW (circles). Right 
panel: summary of the direct and indirect gap for hBN in the different approximations of panels a) and b). 



gies are converged. This procedure yields a quasi-particle 
bandgap for hBN that is about 0.3 eV higher than the 
one obtained on the non-self consistent (Go Wo) level. In 
order to asses how important are the initial eigenfunc- 
tions in this procedure, we performed calculations start- 
ing from LDA or EXX states. 

Electron-hole attraction (excitonic effects) is included 
by solving the Bethe-Salpeter equationii 

{E%-Ell)A s vc ^ Wc , (vck\K eh \v'c'k')A s v , c , k , = tt s A s va 

(0.2) 

are the excitonic eigenstate projections onto the 
electron-hole basis. The interaction kernel K e h "mixes" 
different electron transitions from valence band states 
v,v' to conduction band states c, c' leading to modified 
transition energies fl s . The overall effect in the spectra 
is a redistribution of oscillator strength as well as the 
appearance of new-states within the bandgap (bound ex- 
citons). The GW+BS calculations are performed with 
the code SELB±&. 

In Fig. ^ we present the bandstructure of hBN ob- 
tained using the LDA and EXX approximations and com- 
pare with the GW-calculation. On the LDA and GW lev- 
els, we reproduce quite accurately the results of Ref. U 
The GW-correction increases the direct bandgap (at the 
M-point) by 1.56 eV. Fig. also includes the values for 
the direct and indirect gapsi^. The bare EXX bandstruc- 
ture resembles very closely the LDA GW-bandstructure 
The GW-correction on top of EXX wave- functions yields 
an additional increase of the direct bandgap of 0.45 eV. 
However it still underestimates the "true" direct bandgap 
that must be close to 6.8 eV in order to describe the mea- 
sured optical spectra (see below). 

Our results for the energy-dependent dielectric func- 
tion of hBN (light polarization parallel to the layers) 
are shown in Fig. [2] and compared with the experi- 
mental data of Ref. |6j. The experimental t\ and £2 
were extracted from the electron-energy loss function 
through a Kramers-Kronig analysis 6 . Since the exper- 
imental broadening is about 0.2 eV, we use the same 
value in our calculations. In Fig. [21 a) we present calcu- 



lations based on LDA-wavefunctions. The dash-dotted 
line shows the RPA absorption spectrum (in agreement 
with earlier RPA calculations^). The broad peak with a 
maximum at 5.6 eV is entirely due to the continuum of 
inter-band transitions between the 7r and tt* bands (see 
Fig. [TJ. The calculated GW+BS absorption spectrum 
displays a double peak structure with the main peak at 
5.45 eV and a second peak at 6.15 eV. The nature of the 
spectrum is entirely different from the RPA spectrum: 
the first peak is due to a strongly bound exciton, and 
the second peak contains contributions from higher exci- 
tons and from the onset of the continuum of inter-band 
transitions (see below) . The similarity between the RPA 
and GW+BS spectra stems exclusively from the strong 
broadening employed in the calculation. The main peaks 
in the two spectra are at about the same position because 
of an almost-cancellation between the bandgap widening 
due to the GW-approximation and the red-shift of oscil- 
lator strength due to excitonic effects. Comparison with 
Fig. [21 b) shows that the shape of the GW+BS spectrum 
is in much better agreement with experiment than the 
RPA spectrum. This underlines the importance of ex- 
citonic effects in hBN. The influence of excitonic effects 
becomes even more pronounced when we compare the 
real part of e. Only the GW+BS calculation can repro- 
duce qualitatively the shape of the experimental e±. 

While the shape of the experimental spectrum is well 
reproduced by the LDA GW+BS calculation, the abso- 
lute position is not: The main peak in the experiment 
is 0.65 eV higher than in the LDA GW+BS spectrum. 
We suppose that the quasi-particle gap is higher than the 
one predicted by the LDA GW method. In order to check 
this hypothesis, we performed calculations based on EXX 
wave- functions. In Fig. ^ we showed that the EXX GW 
gap is wider by 0.45 eV than the LDA GW gap. This 
is manifested in Fig. [21 c) where the absorption spectrum 
having similar shape to the LDA GW+BS one is blue- 
shifted by 0.45 eV— . The remaining difference of 0.2 eV 
with respect to the experimental spectrum is most likely 
due to a short-coming of the GW-approximation (pos- 
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FIG. 2: Real (ei) and imaginary (£2) parts of the dielectric 
function of hBN calculated using LDA (a) or EXX (c). The 
experimental data from Ref. is shown in (b) . For the LDA 
results in a) we add the comparison with the results obtained 
at the RPA and full GW+BS levels. The calculations in- 
clude a Lorentzian broadening of 0.2 eV (full-width at half- 
maximum). The light-polarization is parallel to the BN layers. 



FIG. 3: Absorption spectrum of hBN calculated with a 
Lorentzian broadening of 0.025 eV. Solid line: LDA GW+BS 
calculation (with additional scissor operator - see text). Dot- 
ted line: LDA GW RPA. 



sibly due to the neglect of vertex corrections). In what 
follows we will present calculations based on LDA wave- 
function and energies. Beyond the GW-correction, we 
will introduce an additional scissor (0.65 eV) with a value 
such that the first absorption peak at 6.1 eV, in agree- 
ment with experiments^*^ is properly reproduced. 

In Fig.GOwe show the "shifted" GW+BS spectrum with 
a broadening of 0.025 eV. We also display the GW RPA 
spectrum (corresponding to independent quasi-particle 
transitions). The onset of the latter at 6.8 eV marks the 
direct quasi-particle gap with respect to which the ex- 
citon binding energy is calculated. There is a sequence 
of bound excitonic peaks below the onset of the contin- 
uum. Most of the oscillator strength is collected by the 
first bound exciton which is doubly degenerate and has 
a binding energy of 0.7 eV. Upon broadening, the higher 
bound excitons and the absorption at the continuum edge 
form together the second peak at 6.85 eV which has been 
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FIG. 4: Energy scheme of optically active (solid line) and 
"dark" (dashed lines) excitonic states around the first ab- 
sorption peak at 6.1 eV. 



observed in experiments^ and can be seen in Fig. [2] Our 
results are therefore in good agreement with the experi- 
mental data. It is also in good agreement with the recent 
optical absorption measurement in ref. 0, assuming that 
their absorption peak at 6.1 eV is strongly broadened due 
to sample quality of experimental resolution. However, 
the interpretation of the experiments has changed: While 
up to now, the spectrum of hBN has always been inter- 
preted in terms of a continuum of inter-band transitions 
with a direct gap between 5.4 and 5.8 eV, our calculations 
show that the spectrum is dominated by a strongly bound 
discrete Frenkel excitons with a direct quasi-particle gap 
of 6.8 eV. 

Recently, Watanabe et aliiS presented optical absorp- 
tion data that is seemingly in contradiction with previ- 
ous experimental data and with our theoretical results. 
They find evidence for a direct bandgap of 5.971 eV with 
a series of four bound excitons at 5.822, 5.945, 5.962, and 
5.968 eV. The binding energy of the first exciton is there- 
fore 0.149 eV in stark contrast with our value of 0.7 eV. 
However, the absorption data of Ref. [T(J stops at 6 eV, 
making a comparison with previous absorption data diffi- 
cult. We suggest that the peaks that they observe are due 
to the excitation of "dark" excitons, i.e., excitonic states 
that are forbidden by the dipole selection rules but that 
might be accessible when the symmetry of the crystal is 
broken through defects or limited sample quality or when 
spin-orbit interaction is strong enough to lead to spin-flip 
(i.e., excitation of triplet-excitons via intersystem cross- 
ing). In Fig. 01 we present schematically the levels of 
all active and dark excitons in the energy range around 
the first optically allowed exciton at 6.1 eV. At 0.072 eV 
below this value a dark singlet exciton is found, and at 
0.159 and 0.069 eV below the first optically allowed sin- 
glet exciton, two triplet excitons are found. These energy 
values do not correspond exactly to the absorption peaks 
of Ref. [T3, but they may explain qualitatively the occur- 
rence of additional absorption peaks below 6.1 eV in the 
work of ref. 10. The limited accuracy of the approxima- 
tions used in the calculations together with the possible 
role of exciton-phonon coupling could be explain the re- 
maining discrepancies. We note that our calculations are 
at variance with the calculations of Ref. ^| that argue 
that the first absorption peak is composed of 4 different 
optically active excitons. 

Another important check of interpretation is the di- 
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FIG. 5: Calculated loss function of hBN within the RPA (dot- 
ted line) and GW+BS (full- line). Dashed line: experimental 
loss function of ref. 0. Inset: Real and imaginary part the 
dielectric function in the GW+BS and RPA methods. 



rect comparison between calculated and measured loss- 
function of hBN&. This is done in Fig. [3] where we com- 
pare the RPA and GW+BS results. The experimental 
loss function displays a peak around 8.4 eV which is due 
to plasma-oscillations of the 7r-electrons. It has a shoul- 
der at 6.5 eV which has never been explained. A simi- 
lar shoulder has been recently observed in the electron- 
energy loss spectra of multi-wall BN nanotubes^tt. In- 
specting the real part of the calculated dielectric func- 
tion, ex, in the inset of Fig. [21 the origin of this shoulder 



becomes clear. A peak in the loss function occurs close 
to the energies where ex changes its sign from negative 
to positive. While the RPA-calculation yields only one 
such crossing and consequently only one plasmon peak, 
the excitonic effects lead to an additional crossing at 6.5 
eV which causes the shoulder of the 7r-plasmon peak, in 
perfect agreement with experiments. 

In conclusion, we have proposed a new scenario to ex- 
plain the spectroscopic properties of hBN. The main ab- 
sorption peak of hBN at 6.1 eV is due to a strongly bound 
exciton with a binding energy of 0.7 eV. This means that 
the "true" minimum direct quasi-particle gap is 6.8 eV 
which is considerably larger than predicted previously 3 . 
The additional "excitonic" peaks observed in Ref. can 
be explained as due to a relaxation of the selection-rules 
due to symmetry breaking and/or due to spin- flip (sin- 
glet to triplet interconversion) . We also show that the 
electron-energy loss spectra of hBN show a characteris- 
tic influence of excitonic effects that manifest itself as a 
shoulder at lower energies of the main 7r-plasmon peak. 
We expect that these excitonic effects will also play an 
important role in the interpretation of experiments on 
EELS of isolated BN-nanotubesfi. 
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As a consequence of the large EXX-band gap, the static 
dielectric constant, eo = 4.4, calculated in EXX-RPA 
strongly underestimates the experimental value of 4.95—. 
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